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•^T Abstract. We use a simple optical/infrared photometric selection of high 

redshift QSOs which identifies a Lyman Break in the optical and requires a 

red IRAC color to distinguish QSOs from common interlopers. We find 100 

[/-dropout ~ 3) QSO candidates with r' < 22 within 11.2 deg 2 in the ELAIS- 

, ; Nl & ELAIS-N2 fields in the Spitzer Wide-area Infrared Extragalactic (SWIRE) 

^. • Legacy Survey. Spectroscopy of 10 candidates shows that they are all QSOs with 

£f} • 2.83 < z < 3.44. We use detailed simulations which incorporate variations in 

£ — | QSO SEDs, IGM transmission and imaging depth to derive a completeness of 85- 

(f) . 90% between 3.0 < z < 3.4. The resulting luminosity function is two magnitudes 

fainter than SDSS and, when combined with those data, gives a faint end slope 
^D \ /3 = 1.62 ± 0.18, consistent with measurements at z < 2 and steeper than initial 

^vO • measurements at the same redshift. 
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Q ■ 1. Introduction 

a: 

The QSO luminosity function is typically fit to a double power law of the form 
C3 ■ 

>! , v_ 0.92$(M*) 

k^j ; 9(M,Z) - lo0 .4(_ a+1 )( M _ M *) + 10 o.4(-/3+l)(M-A/*) W 

(in magnitudes) where a and (3 are the bright and faint end slopes, respectively. 
At high redshift (z > 2.5), initial indications are that the shape of the QSO 
luminosity function (QLF) may be changing. Hunt et al (2004) (hereafter, H04), 
in a search of faint QSOs at z ~ 3, find a very shallow faint end slope /3 = 
1.24±0.07. At the bright end, Richards et al. (2006, in press) show that the slope 
of the bright end is getting shallower with redshift between 2.5 < z < 5. Given 
the importance of QSO luminosity functions in constrainting galaxy evolution 
models, growth of super-massive black holes (SMBHs), and QSO contribution to 
the ionization of HI and Hell in the IGM, these initial indications of a changing 
QLF warrant larger, deeper surveys to better constrain the high redshift QLF. 
It is difficult to obtain a complete QSO luminosity function at z > 2 for 
two reasons. (1) The colors of QSOs at 2.5 < z < 3.5 coincide with those 
of main sequence stars, resulting in unreliable samples. Therefore, targeted 
follow-up spectroscopy or imaging in redder bands (near-IR) is usually needed to 
distinguish interlopers from QSOs. (2) Surveys of faint QSOs at these redshifts 
require fairly deep imaging (r' ~ 24, for accurate point source determination) 
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over wide areas (> 10 deg 2 ) to find large numbers. The SWIRE Legacy Survey 
(Lonsdale et al. 2003), with IRAC & MIPS (3.6-160/^m) imaging over 49 deg 2 
removes both of these obstacles. 



2. Selection 

The number of Lya absorbers increases rapidly toward higher redshift. The Ly- 
man line and continuum absorption from these systems significantly depresses 
the continuum of high redshift objects blueward of Lya (1216A). Our optical 
color selection looks for this flux decrement (amongst point sources) as it red- 
shifts into the [/-band at z ~ 2.8, causing a red U — g' color in a source with an 
otherwise blue continuum. The color criteria were chosed to select a QSO with 
a 2(j deviation in UV spectral slope and are given in Eq. 1 and plotted in Figure 
1. Typically, high redshift QSO surveys avoid the color-space of main sequence 
stars, but our selection will include many A & F stars as we expect to remove 
these with infrared colors. 

U-g'> 0.33 g'-r'< 1.0 U - g' > 3.9(g' - r') - 2.0 (2) 

As shown by Lacy et al. (2004) and Stern et al. (2005), AGN candidates 
can be selected by their red IRAC colors. Unfortunately, we do not expect 
to detect the majority of r ~ 22 QSOs at these redshifts in the two longest 
wavelength IRAC bands (5.8 & 8.0^m). Therefore, we only use channels 1 & 2 
(hereafter, IRAC1 & IRAC2) since they are two magnitudes more sensitive given 
the same exposure time. In Figure 1 the IRAC colors of all point sources with 
r' < 22 in EN1 are plotted, showing that the IRAC1-IRAC2 color is sufficient 
in distinguishing stars from AGN. Therefore, our infrared criteria are simply 

IRAC1 - IRAC2 > -0.15(AB) /^(4.5/xm) > WfiJy. (3) 

We do not expect this cut to significantly affect our completeness, as the colors 
of our QSO template stay far above this cut between 2.8 < z < 3.5, as do 18 of 
19 SDSS QSOs in the SWIRE survey. 



3. Results 

We search the first two fields for which optical/infrared photometry was avail- 
able (ELAIS-N1 & ELAIS-N2). These fields cover 11.2 deg 2 with both optical 
and infrared imaging, and yield 100 c^-dropout candidates between 19 < r' < 22. 
Ten candidates were followed up with optical spectroscopy with Keck/LRIS and 
Palomar/COSMIC, confirming that all 10 are QSOs with 2.83 < z < 3.44. In 
addition, analysis of the complete optical-IR colors and four band IRAC colors 
(when available) confirm a reliability > 90%. We perform detailed Monte-Carlo 
simulations which account for variations in the QSO spectral slope and emission 
line equivalent widths, the distribution of intervening Lya absorbers, and imag- 
ing depth to derive the selection completeness, C(r' , z), as a function of redshift 
and magnitude. We find that the completeness based on color selection alone is 
85-90% near the center of the redshift window. However, at the faintest magni- 
tudes we are insensitive to some QSOs at the high redshift end (3.3 < z < 3.5) 
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Figure 1. Left panel: The optical color selection. The black solid line is 
the QSO color track and the shaded area is the region occupied when shifting 
the QSO spectral slope by ±2<r. The QSOs are within our color selection 
between 2.9 < z < 3.4. The asterisks are main sequence stars (Gunn & 
Stryker 1983) and the colored dashed lines are color tracks of various GRASIL 
galaxy templates (Silva et al. 1998) between < z < 2 . Right Panel: The 
IRAC colors of point sources (/ < 22) in EN1. The AGN lie at the top right 
of the plot, while the stars occupy the lower left. The AGN selection of Stern 
et al. (2005) (dash-dot) and uur infrared color cut (dashed) are plotted. 



due to the limited depth our [/-band imaging. Therefore, the effective volume 
in the faintest bin is 72% of the volume of the brightest bin. We integrate the 
C(r',z) function to derive an effective volume (and subsequent spece densities) 
for each magnitude bin. 

The binned QLF data are plotted in Figure 2 and match well with the faint 
end values from SDSS (Richards et al. 2006, in press). We fit a double power- 
law to our data and the SDSS data and constrain the bright end slope to the 
value measured by SDSS (a = 2.85). Our fitted parameters are (3 = 1.62 ± 0.18, 
M*(1450A)= -25.6±0.3(AB), and $* = 4.4x 10 7 mag" 1 Mpc" 3 using a ACDM 



cosmology (h = 0.7, Q m = 0.3, ^a = 0.7). 



4. Discussion 



Our fitted faint end slope is significantly steeper than that of H04. A Kolmogorov- 
Smirnov (KS) test gives a probability P=0.004% of measuring a more discrepant 
distribution from the same parent sample. We therefore rule out, at the 4c level, 
the faint slope of H04. We note, however, that at least part of the discrepancy 
may be explained by their use of only broad line (FWHM > 2000km s _1 ) AGN. 
Although we expect most of our objects to have broad lines, we can not rule out 
a significant population of narrow line objects at the faint end of our sample. 
The Pei (1995) fit has an almost identical faint end slope, but with a different 
normalization that predicts a factor of ~ 2 more faint QSOs than our sample. 
This may not be surprising since the shape of the Pei (1995) QLF was deter- 
mined at low redshift and normalized with bright QSOs at high redshift. Since 
the bright end slope gets shallower at high redshift, this causes the faint end of 
the Pei (1995) QLF to be too high. We therefore see no indication of evolution 
of the faint end slope between < z < 3.5. 
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Figure 2. The z ~ 3 QSO luminosity function. The binned data are plot- 
ted from this work (black circles), SDSS from Richards et al. (2006) (blue 
squares), Combo-17 (Wolf et al. 2003) (red diamonds), Hunt et al. (2004) 
(cyan triangles) and Warren et al. (1994) (WHO, green asterisks). Filled sym- 
bols were used in the fitting of the QLF. Also plotted are the fitted QLFs from 
this work (solid), Pei (1995) (dotted), and Hunt et al. (2004) (dashed). 



We have also performed a (/-dropout (z ~ 4) search, with spectroscopic 
follow-up of 10 candidates yielding 7 previously unknown z ~ 4 QSOs. However, 
three sources were galaxies at z ~ 0.4, where the decrement in the g' band 
was a Balmer or 4000A break. Also, this sample suffered somewhwat from 
incompleteness as the redshifting of Ha into IRAC1 at z=3.9 is expected to 
bring many of the QSOs below our IRAC color cut. 

For a more comprehensive analysis of the selection effects, comparison with 
other QLFs and implications for QSO contribution to HI photoionization at 
z = 3, see Siana et al. (2006, submitted). 
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